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Rare semileptonic decays of B and Be mesons are investigated in the framework of 
the QCD-motivated relativistic quark model based on the quasipotential approach. 
Form factors parametrizing the matrix elements of the weak transitions between cor- 
responding meson states are calculated with the complete account of the relativistic 
effects including contributions of intermediate negative energy states and relativistic 
transformations of the meson wave functions. The momentum transfer dependence 
of the form factors is reliably determined in the whole accessible kinematical range. 
On this basis the total and differential branching fractions of the B — t- K^*^l^l^ {uu) 
and Be D^h+r (z/zv), Be — )• D^*^l^l {vv) decays as well as the longitudinal po- 
larization fractions Fl of the final vector meson and the muon forward-backward 
asymmetries ApB are calculated. Good agreement of the obtained results with the 
recent detailed experimental data on the B — )■ K^*^ pL^ ^~ decays from Belle and CDF 
is found. Predictions for the rare semileptonic decays of the Be mesons are given. 

PACS numbers: 13.20.He, 12.39.Ki 

I. INTRODUCTION 

The investigation of the rare weak B and Be meson decays represents a very interesting 
and important problem. Such decays are governed by the flavour- changing neutral currents, 
which are forbidden at tree level in the standard model (SM) and first appear at one- 
loop. Therefore, such decays are very sensitive to the contributions of new intermediate 
particles and interactions, predicted in numerous extensions of the SM (see e.g. [1, 2] and 
references therein). Notwithstanding the fact that such decays have very small branching 
ratios, comparison of existing theoretical and experimental results for the rare semileptonic 
and radiative B — )■ K^*'> decays already provides one of the most rigid constraints on different 
new physics scenarios [3]. 

The theoretical analysis of the rare weak B decays is based on the electroweak effective 
Hamiltonian, which is obtained by integrating out the heavy degrees of freedom (electroweak 
bosons and top quark) [3]. The QCD corrections to these processes due to hard gluon 
exchanges turn out to be important and require to resum large logarithms, which is done 
with the help of renormalization group methods. The operator product expansion allows 
to separate the short-distance part in the B meson decay amplitudes, which is described 
by the Wilson coefficients and can be calculated perturbatively, from the long-distance 
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part contained in the operator matrix elements between initial and final meson states. For 
the investigation of the exclusive decay rates one needs to apply nonperturbative methods 
to calculate these hadronic matrix elements which are usually parametrized in terms of 
covariant form factors. Clearly, such calculation is model dependent. In order to reduce 
model dependence, methods, based on the heavy quark and large energy expansions, have 
been developed. They employ the new symmetries which arise in heavy quark and large 
energy limits and permit to significantly reduce the number of independent form factors 
[4] . Such methods allow a perturbative calculation of QCD corrections to the factorization 
approximation and thus are now popular in the literature [5]. However, the important 
^Qco/i^b corrections cannot be systematically taken into account in such an approach. 

Rare B — )■ K^*"* transitions are the most studied ones both theoretically and experimen- 
tally [3]. Recently, detailed experimental data on differential branching fractions, angular 
distributions and asymmetries in the rare B — > K^*'>iJi^jjr decays became available both 
from B factories and Tevatron [6-9]. The measured values are at present consistent with 
the predictions of the SM within experimental and theoretical uncertainties. Significantly 
better statistics on the rare B decays is expected form LHC experiments (especially from 
LHCb) which will allow precision tests of the SM and can probably reveal signals of new 
physics [10]. It is expected that the Be mesons will be copiously produced at LHC, making 
possible the experimental study of their weak rare decays. Such decays received significantly 
less attention in the literature. The Be D'f^l^l^ and Be D^*^l'^l^ were previously in- 
vestigated using the relativistic constituent quark model [11], light-front quark model [12, 13] 
and three-point QCD sum rules [14] . 

In this paper we study the rare weak B and B,. decays in the framework of the QCD- 
motivated relativistic quark model. Our model was previously successfully applied for the 
investigation of various electroweak properties of heavy and light hadrons. Semileptonic 
decay rates of the B [15] and Be [16, 17] mesons as well as rare radiative decays of the B [18] 
were calculated in agreement with available experimental data. For this purpose, effective 
methods of the calculation of electroweak matrix elements between meson states with a 
consistent account of relativistic effects were developed. They allow to reliably determine 
the form factor dependence on the momentum transfer in the whole accessible kinematical 
range. The form factors are expressed as overlap integrals of the meson wave functions, 
which were obtained in the corresponding calculations of the mass spectra [19, 20]. It is 
important to note that we specially checked [21, 22] the fulfillment of the model-independent 
symmetry relations among form factors arising in the heavy quark and large energy limits. 
Here we apply these methods to the calculation of the form factors of the rare B — > X^*) and 
Be — )■ D^*\D^*^) transitions and on this basis determine branching fractions and differential 
distributions of these decays. 

The paper is organized as follows. The relevant effective weak Hamiltonian for the rare B 
and Be decays is briefly discussed in Sec. II. In Sec. HI we give an outline of our relativistic 
quark model. Then in Sec. IV we discuss the relativistic calculation of the hadronic matrix 
element of the weak current between meson states in the quasipotential approach. Special 
attention is devoted to the contributions of negative energy states and the relativistic trans- 
formation of the wave functions from the rest to the moving reference frame. Form factors 
of the rare semileptonic B K'^*^ and Be D^/^D'^*^) decays are calculated in Sec. V. 
These form factors arc used in Sec. VI for the calculation of the total and differential rare 
decay branching fractions. First we give the necessary formulas and then present our nu- 
merical results. These are then confronted with available experimental data and predictions 
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of other approaches. Finally, Sec. VII contains our conclusions. Expressions for the tensor 
form factors of the rare B and meson decays in terms of the overlap integrals of meson 
wave functions are given in the Appendix. 



II. EFFECTIVE HAMILTONIAN FOR THE RARE B AND Be MESON DECAYS 

The usual approach to the description of rare B decays is based on the low-energy effective 

Hamiltonian, obtained by integrating out the heavy degrees of freedom (the top quark and 
W bosons) of the SM. The operator product expansion separates the short-distance contri- 
butions, which are contained in the Wilson coefficients and can be calculated perturbatively, 
from the long-distance contributions contained in the matrix elements of the local opera- 
tors. The calculation of such matrix elements requires the application of nonperturbative 
methods. 

The effective Hamiltonians for b — )■ fl^l~ and b fvi' transitions {f = s or d), renor- 
malized at a scale /i nib, given by [23] 



i=l 

AC 

= -^v,)VtbCioi, (1) 



where Gp is the Fermi constant, Vtj are Cabibbo-Kobayashi-Maskawa matrix elements, Cj 
are the Wilson coefficients and Oi are the standard model operator basis which can be found 
e.g. in [24]. The most important operators for the b — > fl'^l' transitions are the following 

O7 = ^m,(/a^.(l + 75)6)F'^^ 
O, = ^(/7m(1-75)&)(W, 

Cio = ^(/7m(i - i^mri^i). (2) 

with F^i, being the electromagnetic field strength tensor, and for the b fui' transitions 
we have 

01 = ^(/7m(1 - 75)fe)(^7^(l - 75)^). (3) 
The resulting structure of the free quark decay amplitude has the form: 

M{b ^ fin-) = ^^V,}vJcfih,il - ^.mh^l) + Cio(/7,(l - 75)6)(iV750 



V2 27r 



2m 



g2 



CrUa^,q\l + ^^)b){m) 



M{b ^ fvv) = ^^V,}VtbCl{h,{l - 75)fe)(^7'^(l - 75)^), (4) 

where a is the fine structure constant. 

The effective Wilson coefficient Cf is given [25] by Cf = C7 - C5/3 - Cq, while Cf 
accounts for both perturbative and certain long-distance contributions from the matrix el- 
ements of four-quark operators Ci,...,6- The long-distance (nonperturbative) effects arise 
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from the cc resonance contributions from J / ip ^ ip' . . . and are usually assumed to have a 
phenomenological Breit-Wigner structure. Therefore Cf^ reads as follows [11, 25, 26] 

cf = Cg + J^pert(g') + :^Bw(g')- (5) 

Here the perturbative part is given by 

J^pert(g') = h ^] (3Ci + C2 + 3C3 + C4 + 3C, + Ce) 

\mb ml) 



-\h (^0, (C3 + 3C4) + ^(3C3 + C4 + 3C5 + Ce), (6) 
and the cc resonance part reads 

and is the four-momentum squared of the lepton pair, m^^c are the masses of the b and 
c quarks. The explicit form of the function h{mc/ mb, / m"^) [27] and the values of Wilson 
coefficients Ci,...,io are given in Refs. [11, 25]. 

For the application of the above expressions to the description of the exclusive rare 
semileptonic decays of the B and B^. mesons it is necessary to calculate the matrix elements 
of the operators /7^(1— 75)6 and f(jf^yq^{l—^^)h between initial and final hadron states. Such 
calculation requires the application of nonperturbative approaches. In this paper we use the 
relativistic quark model based on the quasipotential approach for these investigations. 



III. RELATIVISTIC QUARK MODEL 

In the quasipotential approach a meson is described as a bound quark- ant iquark state 
with a wave function satisfying the quasipotential equation of the Schrodinger type 

where the relativistic reduced mass is 

- e;te-, - — UP — ' 

and El, E2 are the center of mass energies on mass shell given by 

— ml + ml — mf + ml , , 
2M ' 2M ■ ^ 

Here M = Ei + E2 is the meson mass, mi,2 are the quark masses, and p is their relative 
momentum. In the center of mass system the relative momentum squared on mass shell 

^2 ^ [M^ - (mi + m2y][M' - (rm - 7712)^] 
^ ^ 4M2 ■ ^ ^ 
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The kernel V{p, q; M) in Eq. (8) is the quasipotential operator of the quark- antiquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by [19] 

V{p, q; M) = u,ip)u2{-p)Vip, q; M)u,{q)u2i-q), (12) 

with 

V(p, q; M) = ^asD,,{kUY2 + Ko„f (k)r;T2;^ + Ko„f (k), 
where as is the QCD coupling constant, Df^^, is the gluon propagator in the Coulomb gauge 

^°°(k) = -g, D^k) = [s^^ - D'^ = D^' = 0, (13) 

and k = p — q. Here 7^ and u{p) are the Dirac matrices and spinors 



where cr and are Pauli matrices and spinors and e{p) — ^/^~+Tn^. The effective long- 
range vector vertex is given by 

r.(k) = 7. + (15) 

where k is the Pauli interaction constant characterizing the long-range anomalous chromo- 
magnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 



e(p) + m 

Hp) 



1 

crp 



e(v) + m 



X 



(14) 



Vv{r) = {l-e){Ar + B), 

Vs{r) = e{Ar + B), (16) 

reproducing 

Konf (r) = Vs{r) + Vvir) ^Ar + B, (17) 

where £ is the mixing coefficient. 

The expression for the quasipotential of the heavy quarkonia, expanded in v'^/c^ can 
be found in Ref. [19]. The quasipotential for the heavy quark interaction with a light 
antiquark without employing the nonrelativistic {v/c) expansion for the light quark is given 
in Ref. [20]. All the parameters of our model hke quark masses, parameters of the hnear 
confining potential A and B, mixing coefficient e and anomalous chromomagnetic quark 
moment k are fixed from the analysis of heavy quarkonium masses and radiative decays. 
The quark masses rrib = 4.88 GeV, rric = 1.55 GeV, = 0.5 GeV, mu,d = 0.33 GeV 
and the parameters of the linear potential A = 0.18 GeV^ and B = —0.30 GeV have the 
values inherent for quark models. The value of the mixing coefficient of vector and scalar 
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FIG. 1: Lowest order vertex function T^^^ contributing to the current matrix element (18). 

confining potentials e — —1 has been determined from the consideration of the heavy quark 
expansion for the semileptonic B D decays [21] and charmonium radiative decays [19]. 
Finally, the universal Pauli interaction constant k = —1 has been fixed from the analysis 
of the fine splitting of heavy quarkonia ^Pj- states [19] and the heavy quark expansion 
for semileptonic decays of heavy mesons [21] and baryons [28]. Note that the long-range 
magnetic contribution to the potential in our model is proportional to (1 + k) and thus 
vanishes for the chosen value of k = — 1 in accordance with the flux tube model. 

IV. MATRIX ELEMENTS OF THE EFFECTIVE WEAK CURRENT 
OPERATORS FOR h ^ s,d TRANSITIONS 

In order to calculate the exclusive rare semileptonic decay rate of the B (Be) meson, it is 
necessary to determine the corresponding hadronic matrix element of the weak operators (2), 
(3) between meson states. In the quasipotential approach, the matrix element of the hadronic 
weak current operator Jj^ , between a B (Be) meson with mass Mb and four- momentum pb 
and a final meson F (F — K^*'> or and D^*'>) with mass Mp and four-momentum pp 
takes the form [29] 

d^pd^q^ 

where r^(p, q) is the two-particle vertex function and ^mpm the meson (M = B,F) 
wave functions projected onto the positive energy states of quarks and boosted to the moving 
reference frame with three-momentum p^- 

The contributions to F come from Figs. 1 and 2. The leading order vertex function 
r^^) corresponds to the impulse approximation, while the vertex function F^^^ accounts for 
contributions of the negative- energy states. Note that the form of the relativistic corrections 
resulting from the vertex function F*^^) is explicitly dependent on the Lorentz structure of the 
quark- ant iquark interaction. In the leading order of the the heavy quark {rub^c oo) and 
large energy expansions for S — )■ F transitions, only F^^^ contributes, while F^^^ contributes 
already at the subleading order. The vertex functions are determined by 

F«(p, q) = Uf{pf)gM(lb){'^^f^{v<i - q.), (19) 

and 

r?'(p,q) = .,(,,)S,W{e.„-j|^,fV(p,-q,) 



(-F(PF)IJf |B(pb)) = / ^^*f-p,(p)r,(p,q)*Bp.(q), (18) 
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FIG. 2: Vertex function F^^^ taking the quark interaction into account. Dashed hnes correspond to 
the effective potential V in (12). Bold lines denote the negative-energy part of the quark propagator. 



+V(p, - q,) 



7i 



lQlt^}ub{(lb)Uq{qq), 



(20) 



where = 7^(1 — 75) for the (axial) vector weak current and = a^^q^il + 75) for the 
(pseudo) tensor current; the subscripts / and q denote the final active s, d and the spectator 
(i, c quarks, respectively; the superscripts "(1)" and "(2)" correspond to Figs. 1 and 2, 
k = p/ - A; k' = qfo + A; A = - p^; 



A(-)(p) = 



e{p) - (m7° + 7°(7p)) 
2e(p) 



Here [29] 



Pb ^ 

qb,q = eb,q{q)^±^n^'\pB)q\ 



i=l 



and n*^*^ are three four-vectors given by 
'«^(p) 



P_ . 



M(£; + M) 



E = Jp^ + M2. 



It is important to note that the wave functions entering the weak current matrix element 

(18) are not in the rest frame in general. For example, in the B meson rest frame (p^ = 0), 
the final meson is moving with the recoil momentum A. The wave function of the moving 
meson "^p/^ is connected with the wave function in the rest frame "^fo = by the 
transformation [29] 

*^a(p) = Df(RYjDl/\RYj^po{p), (21) 

where is the Wigner rotation. La is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D^/'^[R) in spinor representation is given by 



(22) 
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where 



'{p) + m 



1 + 



ap 



2m \ ' e(p) + 
is the usual Lorentz transformation matrix of the four-spinor. 



V. FORM FACTORS OF RARE SEMILEPTONIC DECAYS 

The matrix elements of the weak current for rare B decays {B denotes either B or B^.) to 
pseudoscalar mesons (P = K,Ds, D) can be parametrized by three invariant form factors, 



„ Ml -Ml , ■ 



(23) 



{P{pF)\q<j^^qMB{pB)) 



Mb + Mp 



[q'ip'B+P'F)-{Ml-Ml)q% 



(24) 



where /+(0) = /o(0), g = Pb—Pf, and M^^p are the masses of the B meson and pseudoscalar 
P meson, respectively. 

The corresponding matrix elements for the rare B decays to vector mesons {V — 
K*,D*, D*) are parametrized by seven form factors. 



{V{pF)\qrh\B{pB)) 
{V{pF)\qri,h\B{pB)) 



2iV{q^) 
Mb + Mv 



e'"'"'e:pBpPFa: 



(25) 



q 



2MvAo{q') — 

q 

e -q 



*fj. 



-A2{q' 



qi^ + {MB + Mv)A,{q^){e 



Ml- M?, ,' 
P^b + Pf- ^ , q^ 



^q^ 



Mb + Mv 

{V{pF)\qia^-q,h\B{pB)) = 2T^{q^)e^-P'^ elpppPB., 
{V{pF)\qia^^^,qMB{pB)) = Uq^)[{Ml - Ml)e*^ - {e* ■ q){p% + p%)\ 



(26) 
(27) 



+Uq'W-q) 



Ml 



Ml 



{Pb+Pf) 



(28) 



where 2MvAo{0) = {Mb + My)Ai(O) - {Mb - Mv)A2{0), Ti{0) = r2(0); My and are the 
mass and polarization vector of the final vector meson. 

We previously studied the form factors {f+,fo,V,Ao,Ai,A2) parametrizing the matrix ele- 
ments of vector and axial vector charged weak currents for B — )■ 7r(p) [15] and Be — )■ r]c{J/i^), 
Be — )■ D^*"* [16], B(. — )> B^*\B^*^) [17] transitions in the framework of our model. The neces- 
sary formulas for these form factors can be found in Appendix of Ref. [16]. Now we apply 
them to the calculation of the form factors, parametrizing neutral current matrix elements 
for the B — > K'^*^ and Be — > D^*\ Be — > D'^*^ transitions. For the remaining tensor form 
factors we use the same approach described in detail in Rcfs. [15-17]. Namely, we calculate 
exactly the contribution of the leading vertex function T^^^ (19) to the transition matrix 
element of the weak current (18) using the 5-function. For the evaluation of the subleading 
contribution F^^^ we use expansions in inverse powers of the heavy 6-quark mass from the 
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TABLE I: Form factors of the rare semileptonic decays B — )■ K^*^l'^l calculated in our model. 
Form factors friq^)-, Viq^)-, ^0(9^), Ti{q^) are fitted by Eq. (29), and form factors foiq^), 

^i(g2), A2{q''), T2{q% n{q^) are fitted by Eq. (30). 
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initial B meson and of the large recoil energy of the final heavy-light meson. Note that the 
latter contributions turn out to be rather small numerically. Therefore we obtain reliable 
expressions for the form factors in the whole accessible kinimatical range. It is important 
to emphasize that doing these calculations we consistently take into account all relativistic 
corrections including boosts of the meson wave functions from the rest frame to the moving 
one, given by Eq. (21). The obtained expressions for the tensor form factors /j',Ti,T2,T3 
are presented in Appendix (to simplify these expressions the long-range anomalous chromo- 
magnetic quark moment was explicitly set as k = —1). In the limits of the infinitely heavy 
quark mass and large energy of the final meson, the form factors in our model satisfy all 
model independent symmetry relations [4, 22]. 

For numerical calculations of the form factors we use the quasipotential wave functions 
of the B, Be, K*, Dg and D mesons obtained in their mass spectra calculations [19, 20]. 
Our results for the masses of these mesons are in good agreement with experimental data 
[30], which we use in our calculations. 

We find that the rare semileptonic B K^*H'^l^ and Be — )■ D^*\D^*^)l^l~ decay form 
factors can be approximated with good accuracy by the following expressions [15, 31]: 

(a) F(e) = {U(l\h(d\V(q\A„(e),T,(e)} 

= 7 / "^'"'^^ (29) 

(b) F(<(=) = [h(q'),A,{q^),A^(q%T,{q^),T,(q-')} 

where M = Mb^ for Aq and M = Mb* for all other form factors (for Be — )■ D^*H~^l~ decays 
M (*) should be replaced by Mb(*))- The values F(0) and ui 2 are given in Tables I-III. The 
difference of fitted form factors from the calculated ones does not exceed 1%. We plot these 
form factors in Figs. 3 and 4. 




VI. RESULTS AND DISCUSSION 

Now we use the obtained form factors for the numerical calculation of decay rates and 
other important observables of the rare semileptonic B decays and confront their values with 
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TABLE II: Form factors of the rare semileptonic decays Be — s- Ds l^l calculated in our model. 
Form factors friq^), Vio^), Ti{q^) are fitted by Eq. (29), and form factors fo{q'^), 

Ai{q^), A2iq^), T2{q^), T^{q^) are fitted by Eq. (30). 
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TABLE III: Form factors of the rare semileptonic decays B^ — s- D^*H^l calculated in our model. 
Form factors f+{q^), friq'^), Vig'^), ^o('7^), Ti{q'^) are fitted by Eq. (29), and form factors /o(g^), 
^i((?2), A2iq^), T2{q^), niq^) are fitted by Eq. (30). 
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FIG. 3: Form factors of the B — t- K^*^ decays. 
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FIG. 4: Form factors of the Be ^ Ds decays. 



available experimental data. 



A. B ^ Pl+l- and B Vl+l' decays 



The matrix element of the b — )■ fl^l (/ = s or d) decay amplitude (4) between meson 
states can be written [11, 25] in the following form 



M{B pi+r 

M(B Vl+l 



Gpa 
2^271 



et'^r«(/7'^/) + et'^T£)(/7^5/) 



(31) 



where T*-*) are expressed through the form factors and the Wilson coefficients. Then these 
amplitudes can be written in the helicity basis e^{m) as (see [11]) 
(a) B ^ P transition: 

i/«=etM(^)TW, (32) 

where 



H. 



± 
(1) 



0, 

Ai/2 



2mh 



Mb + Mp 



12 



(2) 



(2) 



AV2 
Ml 



Ml - Ml 



(33) 



Here A = A(M|, Ml, q^) = M^ + Mj, + - 2{MlMl + Mlq^ + Mlq^) and the subscripts 
±, 0, t denote transverse, longitudinal and time helicity components, respectively, 
(b) S — >■ y transition: 

i/»=et/^(m)et^T», 



(34) 



where e'^ is the polarization vector of the vector V meson and 



^(2) 



(1) 



(2) 



-{Ml - M'y) 



± — 



a 



eff 



V{q^) 



Mb - My 
2mb 



+ ^-^C^/^Uq^) 



Mb + Mv 



+ 



C^/%{q') 



iMB^MM)±X^^^j£^^ 



2My^/? 
1 



CI" |(Ml - - q^){MB + My)A,{q') - ^^^Mq") 



{Ml + 3M'y - q')T,{q') 



A 



Ml 



Ml 



Uq') 



c 
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2MyvV 
^"\cl"A,{q% 



{Ml -M^- q^){MB + Mv)A^{q^) 



A 



Me + My 



Ai/2 



Cio^o(g'), 



(35) 



here A = A(M|, , g^) = m| + + - 2(M|M2 + M^g^ + 
The differential decay rate then reads [11] 



dT{B P{V)l+l-) 



2% 



48M|\ 



+^(2)^t(2) ( ] 



+ 



4mf ' 



(36) 



where m; is the lepton mass and 



'0 -'-'0 



(37) 



The forward-backward asymmetry is given by 

Re{H^^^Hf^)-Re{H^}^H^}'^) 



^FB 



3 

4^ 



^ //(l)//t(l) f 1 + l!^ + i/(2)//t(2) f 1 _ 



(38) 
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The longitudinal fraction of the V polarization has the form 



o-{l)o-t(l) 



1 + ^1+ H^^'H^ ( 1 - ^ ) + "^m^'Hi 



r(2) o-t(2) 



Ami 



2mf 



r(2) o-t(2) 



HWHm ( 1 + ^ ) + //(2)//t(2) f 1 - i!# ] + '-^sh^'^h}^'^ 



(39) 



These two observables are the most popular quantities for B K*{^ K7T)ii^fi^ decays, 
since they are convenient for the experimental measurements. They enter the decay differ- 
ential distributions in cos 9k 



IdViB ^ K*a+a-) 3^ 3, o.. 



(40) 



and in cos On 



1 dr(S K*fx+fj,-) 3 



d cos 



= -Fi(l-cos^^^) + -(l-F^)(l + cos'^^) + Ai.BCOs^^, (41) 



where 9k is the angle between the kaon direction and the direction opposite to the B meson 
in the K* rest frame, and 9^ is the angle between the and the opposite of the B direction 
in the dilepton rest frame. Therefore they can be determined experimentally using the 
angular analysis. 



B. B ^ Pvu and B Vuv decays 
The differential decay rate for the B — > P{y)vv is given by 



= 3- 



tf vtb\ 



X 1/2^2 

24M| 



dq^ (27r)3 V 27r 

where the factor of 3 originates from the sum over neutrino flavours, 

and the helicty amplitudes H!^^ read as follows 
(a) B ^ P transition: 



± 



0, 



(h) B ^ V transition: 



H 



-{Mb + My)A,{q') ± X^^^^^^l^ 



2MvVq 



X 



(42) 



(43) 



Mb + Mx 



-A,{q' 



(44) 
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TABLE IV: Comparison of our predictions for the rare semileptonic B — >• K^*^ decay nonresonant 
branching fractions with experimental data (in 10~^). 

Decay our BaBar [6, 33] Belle [7, 34] CDF [8] CDF [9] HFAG [35] 

B K^+^- 4.19 3.4 ±0.7 ±0.2 4.8^^^ ± 0.3 5.9 ± 1.5 ± 0.4 3.8 ±0.5 ±0.3 4.5 ±0.4 
B Kt+t- 1.17 

B Ku9 26.1 < 140 

B^K*ii+fi- 9.25 7.8li;?±l.l 10.7li;J ± 0.9 8.1 ±3.0 ±1.0 10.6 ±1.4 ±0.9 10.8li:i 

B K*T+T~ 1.03 

B K*vv 63.2 < 800 



Here Q = -X{xt)/ sin^ Ow, with Xt = m^/m'^, 9\Y is the Weinberg angle, and the function 
X{xt) at the leading-order in QCD has the form 

, , X [2 + X 3x - 6 , \ 

while the next-to- leading order expressions are given in Ref . [32] . 

Substituting the current experimental values for the top {rrit) and VF-boson (rnw) masses 
one gets [2] 

CI = -6.38 ± 0.06, (45) 

where the error is dominated by the top quark mass uncertainty. In the following calculations 
we use the central value of C^. 

The differential longitudinal polarization fraction Fl of the V meson is defined similar to 
Eq. (39) 



(46) 



C. Numerical results 



Now we substitute the above calculated form factors in the expressions for decay rates, 
asymmetries and polarization fractions and perform numerical calculations. 

First we compare the predictions of our model for the B — )■ Kl^l~ and B — )■ K*l^l~ 
decays with available experimental data. The calculated values for the branching fractions 
of the rare semileptonic decays B — )■ K^*H'^l~ and B — )■ K'^*^i'V and available experimental 
data are given in Table IV. We find good agreement of our results for the B — > Kji^ji' 
and B — > K*fi~^fi~ decays with experimental data. A more stringent test of our predictions 
can be achieved by comparison with new data on differential decay distributions. This is 
done in Fig. 5 where we confront our predictions for differential decay rates, the longitudinal 
polarization fraction of the K* and the muon forward-backward asymmetry Aps with 
detailed experimental data from Belle [7] and CDF [9]. In this figure wc plot our results 
both for the nonresonant (solid line) quantities and quantities including the J/i/j and i// 
resonance contributions (dashed line). Note that the resonance regions are vetoed in the 
experimental analysis. Reasonable agreement of our predictions with experimental data is 
found. The current experimental data on ApB are not precise enough to give a definite 
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FIG. 5: Comparison of theoretical predictions for the differential branching fractions dBr/dq^, 
the K* longitudinal polarization Fl and muon forward-backward asymmetry Aps for B — )• K^*'^ 
decays with available experimental data. Nonresonant and resonant results are plotted by solid 
and dashed lines, respectively. Belle data are given by dots with solid error bars, while CDF data 
are presented by filled circles with dashed error bars. 



conclusion whether this asymmetry has a zero or not. Our model predicts the zero oi ApB 
at ql = 2.74 GeV^ which is in agreement with the value ql = 2.88lo;28 GeV^ given in [25]. 
It is expected that the accuracy of experimental data will increase significantly in the near 
future. 

In Fig. 6 we plot our results for the differential branching fractions of the B — > Kt^t~ 
and B — )■ K*t^t^ decays. The calculated values for these decay branching fractions are 
presented in Table IV. There we also give our results for the B — )■ Kuu and B — )■ K*!/^ 
branching fractions. None of these modes have been measured yet. Only experimental 
upper bounds have been recently set on branching fractions for the B — > K*^!) decay by 
BaBar [33] and for the B — )■ Kvv decay by Belle [34]. These bounds are about an order 
of magnitude higher than our model predictions. In Fig. 7 we show our predictions for 
the differential branching fraction and the K* longitudinal polarization fraction Fl for the 
B — )■ K*^ decay. As it is noted in Ref. [2], the value -Fl(O) = 1 is imposed by helicity 
conservation, while Fi{q^^-^ = 1/3 follows from the absence of a preferential direction 
at the point q^s.x — {^b — Mk*Y where both K* and B are at rest. The differential 
branching fraction of the B — )• K*i>u decay as well as its integrated value are close to the 
ones {Br{B — )■ K*uu) = (6.8li;?) x 10~^) given in Ref. [2], while the shape of Fl is slightly 
different. We also get the value for the q^ integrated longitudinal polarization fraction 
(^Fl) ^ 0.31, which is significantly lower than the one of Ref. [2] {F^) = 0.54 ± 0.01. Note 
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FIG. 6: Predictions for the differential branching fractions of S — )• K^*^t^t decays. Nonresonant 
and resonant results are plotted by solid and dashed lines, respectively. 
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FIG. 7: Predictions for the differential branching fractions and the K* longitudinal polarization 
fraction Fl for the B — >■ K*vu decay. 



that our results for the branching fractions of the B — )■ K^*'>vi' decay, though consistent 
with the ones from [2], are somewhat lower than the predictions from [23, 26]. 

Next we present our results for the rare semileptonic Be decays. In Table V we com- 
pare available theoretical predictions for the branching fractions of the rare semileptonic 
Be — )■ D^*^ and Be — )■ D^*^ decays. The investigations [11-13] are based on the relativistic 
constituent quark model and light-front quark models, while the authors of Ref. [14] use 
three-point QCD sum rules. Although the results of these approaches are consistent in the 
order of magnitude of the branching fractions, they differ by more than a factor of 2 for 
some decay modes. We find the best overall agreement of our predictions for the branching 
fractions with the results of the relativistic quark model [11]. The differential branching 
fractions, the longitudinal polarization Fl of the final vector meson and the muon forward- 
backward asymmetry ApB for the B^ — )■ Dg and Be — )• -D* transitions in our model are 
plotted in Figs. 8 and 9. Similar curves have been obtained for Be — )■ -D*^*^ transitions 
and are not shown here. We predict the following values of the position of the zero of 
the forward-backward asymmetry ApB'- Qq = 2.41 GeV^ for the Be — DljX^ ji' decay and 
ql = 2.46 GeV2 for the Be D*fi+fi- decay. 
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TABLE V: Comparison of theoretical predictions for the nonresonant branching fractions of the 



rare 


semileptonic Be - 


-> Dl and Be - 


> decays (m 
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Decay 


our 


[11] 


[12] 


[14] 


[13] 


B - 


-4- D //+;y- 


11.6 
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6.1 ± 1.5 


5.4 




D.T^T^ 
~ ^ s ' ' 
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65 
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FIG. 8: Predictions for the differential Be — s- Ds decay branching fractions. Nonresonant and 
resonant results are plotted by solid and dashed lines, respectively. 



VII. CONCLUSIONS 

In this paper we obtained the form factors of rare semileptonic decays of the B and 
Be mesons in the framework of the QCD-motivated relativistic quark model based on the 
quasipotential approach. The consideration is done with a systematic account of all rel- 
ativistic effects, which are very important for such transitions. Particular attention was 
devoted to the inclusion of negative-energy contributions and to the relativistic transforma- 
tion of the meson wave function from the rest to the moving reference frame. As a result, 
the dependence of these form factors was explicitly determined in the whole accessible 
kinematical range without using any ad hoc assumptions and extrapolations. It is important 
to point out that the obtained form factors satisfy all heavy quark and large energy sym- 
metry relations in the corresponding limits [22]. Note that the resulting decay form factors 
are expressed through the overlap integrals of the initial and final meson wave functions. 
The relativistic wave functions, obtained previously in the investigations of the meson mass 
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FIG. 9: Predictions for the differential brandling fractions dBr/dcf', the longitudinal polarization 
Fl and muon forward-backward asymmetry ApB for Be — )• decays. Nonresonant and resonant 
results are plotted by solid and dashed lines, respectively. 



spectra, were used for the numerical calculations. This significantly improves the reliability 
of the calculated form factors. On the basis of these form factors branching fractions and 
different differential decay distributions were obtained. 

First we tested our model by confronting its results for the B — Kfi^n~ and B 
K*^'^^~ decays with the available detailed experimental data. It was found that the total 
and differential branching fractions, the K* meson longitudinal polarization fraction and 
the muon forward-backward asymmetry Aps agree well with data. 

Secondly, we presented detailed predictions for the rare semileptonic decays of the Bf. 
meson which can be investigated in the LHCb experiment at CERN, where the Be mesons 
are expected to be copiously produced. Finally, we compare our results on these decays 
with the ones previously available in the literature in Table V. The predictions for the 
differential branching fractions, the vector meson longitudinal polarization fraction Fl and 
muon forward-backward asymmetry ApB are also given in Figs. 8 and 9. 
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Appendix: Tensor form factors of rare B and Be decays 

(a) B P (B K, Be Dg and B^-^ D) transiUons (see Eq. (24)) 
Mq') = A'\q') + ef?'\q') + (1 - e)ff^'\q% 



(A.l) 



X • 
2 

+- 



e/(p + A) +m/ 



+ 



(pA) 

\^e/(p + A) + m/ eb(p) 
1 1 



+ mbj 



X 



3 + Mp + A) + m/ 

/I 1 

+ 



^e/(p + A)+m/ e6(p) + m6 



U*b(p), 



(A.2) 



lEp 
Mb 
1 



Ep + Mp )\ 2e/(p + A) 



, 6/(p + A)-mf (pA) \ 

2e/(p + A)(e/(p + A) +mj) V 2m6 A^^ 



X 



Mp - e/ p + 



2e„ 



Ep + Mp J ' Ep + Mp 



2e„ 



(pA) [ Mb + - 6,(p) - 6,(p) - 6^ (p + ^^A) - 6, (p + ^^ A) 
A2 [ 2e/(p + A)(e/(p + A)+m/) 

^ Mb - Mp - 6,(p) - 6,(p) + 6^ (p + a) + e, (p + a) 



2mb(eb(p + A) + m^) 



■^b(p), 

(A.3) 



= (Mb + Mp)i 

(pA) r 



P 



2e„ 



X- 



Mp y (27r)3 ^ " V ' ^P + Mp ; \| 2e/(p + A) 
Mp - eb(p) - eq{p) 



e/(p + A) +m/ 



A2 \2ey(p + A)(e^(p + A)+m^) 

1 / 1 e/(p + A)-m/ 



2m6 Ve6(p + A) + mb 2e/(p + A)(e/(p + A) + m/) 



X 



Mp - e/ p + 



2e„ 



£;p + Mp 



2e. 



Ep + Mp 



>*p(p), (A.4) 



where the superscripts "(1)" and "(2)" correspond to Figs. 1 and 2, £ is the mixing coefRcient 
in the confinement potential indicated in Eq. (16) and 



|A| = 



{Ml + M| - ?2)2 



4M| 



M|, 
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Ep = -\/m|, + A2, eqip + aA) = + (p + aAf {Q ^ b,c,s,u,d). 

Here S stands for the B or meson, P — K, Dg, D is the final pseudoscalar meson, f — s,d 
is the final active quark and q — u,d,c denotes the corresponding spectator quark. 

(h) B (B ^ K*, B^-^ D* and B^ -> D*) transitions (see Eqs. (27) and (28)) 

TM') - Ti'\q^) + sT,'^'\q^) + (1 - e)Tr^'\q% (A.5) 



Ti'\q') 



E- 



V 



x|l + 



d^p 



P + 



2e„ 



(27r)3 ^ " V" ' + My )\ 2ef{p + A) \ 2e,{p) 
pV3 + (pA) 



ef{p + A) + 



65 (p) + TTlfe 



(e/(p + A) + mf){eb{p) + rrib) 
p2A2 



3(£;y + My)(eg(p) + mg)(e/(p + A) + mf){ek{p) + m,,) 



+ (Mb - Ev] 

^2 



€,(;> + A) 
(- ' 



+ 



1 



+ 



1 



(pA) / 

A2 i^e_^(p + A) +m/ ' eb{p) + rUb J 



+ 



3{Ev + My) \eq{p) + ruq \ebip) + e/(p + A) + mf J 

2 



(e/(p) +m/)2 



^*b(p), 



(A.6) 



d^p - 
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2e„ 



e/(p + A) -ruf 



e/(p + A) +m/ 
^ 2e;(p + A) 
Mb - Ev (pA) ■ 



\2ef{p + A){ef{p + A)+mf)\ 2mb A?) 

X (^^ - ray^) - ^^)) - - ^^^^ 

/ Mb + My - 6,(p) - 6,(p) - 6/ (p + ^^A) - 6. (p + ^^A) 
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(A.7) 
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X 



2m6(eb(p + A) + m^) 



*s(p), (A.8) 



2e/(p + A)(e/(p + A)+m/), 

T2(g^) = ri^^g^) + erf + (1 - e)T^^'\q% (A.9) 
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\ 2m,(e,(p + A)+m,) 

Mb - eb{p) - eg{p) 



^b(p), (A.12) 



2e/(p + A)(e/(p + A)+m/)^ 
Uq') = T^iq') + eT,'^'\q') + (1 - e)Tp\q'), (A.13) 
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